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We consider the reactions 7p -^ MB, with MB being either K+ A, K+ E", or 
tt"*" n, within a diquark model which is based on perturbative QCD. The model pa- 
rameters and the quark-diquark distribution amplitudes of the baryons are taken 
f-f^ ' from previous investigations of electromagnetic baryon form factors and Compton- 

scattering off protons. Reasonable agreement with the few existing photoproduc- 
tion data at large momentum transfer is found for meson distribution amplitudes 
compatible with the asymptotic one (oc x{l — x)). We present also first results for 
hard electroproduction of the K+ A final state. Our predictions exhibit some char- 
acteristic features which could easily be tested in forthcoming electroproduction 
experiments. 
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0^ ' Our investigation of exclusive photo- and electronmduction of mesons is part 

of a systematic study of hard exclusive reactions crBa within a model which is 
Q_i' based on perturbative QCD, in which baryons, however, are treated as quark- 

' , diquark systems. This model has already been applied to baryon form factors 

g^' in the space- El and time-like region El, real and sdrtual Compton scattering 0, 

r^ , two-photon annihilation into proton-antiprotona and the charmonium decay 

77c —> p pl3. A consistent description of these reactions has been achieved in the 
sense that the corresponding large momentum-transfer data (p^>,3 GeV^) are 
p\^ ' reproduced with the same set of model parameters. Further applications of this 

model include three-body J/^ decays Q and the calculation of Landshoff con- 
tributions in elastic proton-proton scatteringQ. Like the usual hard-scattering 
approach (HSA)la the diquark- model relies on factorization of short- and long- 
distance dynamics; a hadronic amplitude is expressed as a convolution of a 
hard-scattering amplitude T, calculable within perturbative QCD, with distri- 
bution amplitudes (DAs) (j)^ which contain the (non-perturbative) bound-state 
dynamics of the hadronic constituents. The introduction of diquarks does not 
only simplify computations, it is rather motivated by the requirement to ex- 
tend the HSA from (asymptotically) large down to intermediate momentum 
transfers (p^>,3 GeV^). This is the momentum-transfer region where some 
experimental data already exist, but where still persisting non-perturbative 
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effects, in particular strong correlations in baryon wave functions, prevent the 
pure quark HSA to become fully operational. Diquarks may be considered as 
an effective way to cope with such effects. Photoproduction of mesons is one 
class of exclusive reactions wher£ the pure-quark HSA obviously fails in repro- 
ducing the hard-scattering dataQ and where a modification of the perturbative 
QCD approach seems to be necessary in the few-GeV region. 

The following section starts with an outline of the hard-scattering approach 
with diquarks and introduces the hadron DAs to be used in the sequel. The 
diquark model predictions for the various reaction channels are presented in 
Sec. H along with some general considerations on photo- and electroproduction 
reactions. Our conclusions are finally given in Sec. m 

2 Hard Scattering with Diquarks 

Within the hard-scattering approach a helicity amplitude Mi^} for the reaction 
7^*' p ^ MB is (to leading order in 1/p_l) given by the convolution integrals 

M{\}{s,i)= dxidyidzKJ)^ (zi,pj_)0^ (yi,pj.)f{A}(xi, j/i, zi; s, t)0P(a;i,p_i) . 
Jo 

(1) 
The distribution amplitudes (ff^ are probability amplitudes for finding the va- 
lence Fock state in the hadron H with the constituents carrying certain fractions 
of the momentum of their parent hadron and being coUinear up to a factor- 
ization scale p±^ . In our model the valence Fock state of an ordinary baryon is 
assumed to consist of a quark (q) and a diquark (D). We fix our notation in 
such a way that the momentum fraction appearing in the argument of 0^ is 
carried by the quark - with the momentum fraction of the other constituent 
(either diquark or antiquark) it sums up to 1 (cf. Fig. nl). For our actual cal- 
culations the (logarithmic) p±^ dependence of the DAs is neglected throughout 
since it is of minor importance in the restricted energy range we are interested 
in. The hard scattering amplitude TfA} is calculated perturbatively in coUinear 
approximation and consists in our particular case of all possible tree diagrams 
contributing to the elementary scattering process 7qD -^ qqqD. A few exam- 
ples of such diagrams are depicted in Fig. H The subscript {A} represents the 
set of possible photon, proton and A helicities. The Mandelstam variables s 
and t are written with a hat to indicate that they are defined for vanishing 
hadron masses. The calculation of the hard scattering amplitude involves an 
expansion in powers of (?tih/s) which is cut off after ©(ttih/s). Hadron masses, 
however, are fully taken into account in fiux and phase-space factors. 

The model, as applied in Refs. 1-3, comprises scalar (S) as well as axial- 




Figure 1: Graphical representation of the hard-scattering formula, Eq. (hh, for 7p — > MB. 
Xi, yi, and Zi denote longitudinal momentum fractions of the constituents, the Ai's their 

helicities. 



vector (V) diquarks. V diquarks are important, if one wants to describe spin 
observables which require the flip of baryonic hehcities. The Feynman rules for 
electromagnetically interacting diquarks are just those of standard quantum 
electrodynamics. Feynman rules for strongly interacting diquarks are obtained 
by replacing the electric charge e by the strong coupling constant gs times 
the Gell-Mann colour matrix t". The composite nature of diquarks is taken 
into account by multiplying each of the Feynman diagrams with diquark form 
factors F^ (Q^) which depend on the kind of the diquark (D = S, V), the 
number n of gauge bosons coupling to the diquark, and the square of the 4- 
momentum Q^ transferred to the diquark. The form factors are parameterized 
by multipole functions with the power chosen in such a way that in the limit 
Pj^ -^ oo the scaling behaviour of the pure quark HSA is recovered. 

Having outlined how diquarks are treated perturbatively, a few words 
about the choice of the quark-diquark DAs (which incorporate the q-D bound- 
state dynamics) are still in order. In Refs. 1-3 a quark-diquark DA of the form 

(ci — C2 — for S diquarks, x longitudinal momentum fraction carried by 

the quark) 



(j)-£,{x) oc x{l — x) (1 + cix + C2X ) exp 



-6 



2 / "^q 



[l-x] 



, D = S,V, 

(2) 

in connection with an SU(6)-like spin-flavour wave function, turned out to be 
quite appropriate for octet baryons B. The origin of the DA, Eq. (|^), is a nonrel- 
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Figure 2: A few representative examples for three-, four-, and five-point contributions to 

7P ^ MB. As indicated, diagrams are first calculated in Born order for point-like diquarks 

and afterwards multiplied with appropriate vertex functions (diquark form factors). 



ativistic harmonic-oscillator wave functions. Therefore the masses appearing 
in the exponential have to be considered as constituent masses (330 MeV for 
light quarks, 580 MeV for light diquarks, strange quarks are 150 MeV heavier 
than light quarks). The oscillator parameter b^ = 0.248 GeV~^ is chosen in 
such a way that the full wave function gives rise to a value of 600 MeV for the 
mean intrinsic transverse momentum of a quark inside a nucleon. 

For meson DAs various models can be found in the literature. In order 
to check the sensitivity of our calculation on the choice of the meson DAs we 
employ two qualitatively different forms. On the one hand, the asymptotic DA 

(/iasy(a;) ex a;(l -x) , (3) 

which solves the p±_ evolution equation for (j)(x,pj_) in the limit p± — > oo, and, 
on the other hand, the DAs 

<^5z(a;) ex <^asy (2a; - 1)" , 0gz(a;) « <^asy [0-08 + 0.6(2x - l)^ - 0.25(2.x - l)^] , 

(4) 
which have been proposed in Ref. 9 on the basis of QCD sum rules. The 
"normalization" of the meson DAs is determined by the experimental decay 
constants for the weak vr, K ^ ^v^^ decays. The analogous constants /s and 
/v for the q-D DAs of baryons are free parameters of the model. They are, in 
principle, determined by the probability to find the q-D state (D = S, V) in the 
baryon B and by the transverse-momentum dependence of the corresponding 
wave function. 

For further, details of the diquark model we refer to the publication of 
R. Jakob et al.lll. The numerical values of the model parameters for the present 
study are also taken from this paper. 



3 Photo- and Electroproduction Reactions 

Exclusive photoproduction of pseudoscalar mesons can_iii general, be described 
by four independent helicity amplitudes M\^^XB,X-yAp^ 

D = Mq +i +1 _i , 52 = Afg^+i _^i +1 . (5) 

N , Si, S2, and D represent non-flip, single- flip, and double- flip amplitudes, 
respectively. Two additional amplitudes, 

N = il/o_+i^o,+ i , 53 - Mo,_i,o,+i > (6) 

occur if electroproduction is considered. Electroproduction and photoproduc- 
tion are related as far as in the one-photon approximation the electroproduc- 
tion cross section can be expressed in terms of four response functions dcru? 
daL, (icTx, and dai fpr the production of a pseudoscalar meson by a virtual 
photon, 7* p ^- MBO. The two helicity amplitudes in Eq. (0) correspond to 
the longitudinal polarization degree of the virtual photon. The two response 
functions da\j and duL are just cross sections for (spin averaged) transversely 
and longitudinally polarized photons, respectively. The other two functions 
are transverse-transverse (dcT) and longitudinal-transverse {dai) interference 
terms. In the limit of vanishing photon-virtuality {Q^ -^ 0) dan reduces to 
the ordinary photoproduction cp^ss section, the ratio daT/dau becomes the 
(negative) photon asymmetry EEj and dai, as well as dai become zero. 
From the q-D flavour functions (D = S, V) 

Xs=uS[u,d], Xv = [uV{u,d} - \/2dV{u,u}]/V3, (7) 

X§ = dS[u,d], X^ = -[dV{,,d}-%/2uV{d,d}]/V3, (8) 

xf = [dS[„,,] +uS[d^s]]/V2, xv =- [2sV{,,d} - dV{,,3} - uV{d,s}]/%/6 , (9) 

Xs = [uS[d,s] - dS[u,s] - 2sS[u,d]]/V6 , Xv = [uV{d,s} - dV^u^s}]/V2 ■ (10) 

for the baryons one infers already that the three production channels we are 
interested in differ qualitatively in the sense that K+ A is solely produced 
via the Sfu d] diquark, K+ E" via the V/u d} diquark, and 7r+ n via S and 
V diquarks. An important consequence of this observation is that helicity 
amplitudes and hence spin observables which require the flip of the baryonic 
helicity are predicted to vanish for the K+-A final state (e.g., the polarization 
P of the outgoing A). For the other two channels helicity flips may, of course, 
take place by means of the V diquark. 



One of the qualitative features of the HSA is the fixed-angle scaling be- 
haviour of cross sections. The pure quark HSA implies an s~^ decay of the 
photoproduction cross section. This scaling behaviour is, of course, recovered 
within the diquark model in the limit s — > cx). However, at finite s, where the 
diquark form factors become operational and diquarks appear as nearly ele- 
mentary particles, the s~^ power-law is modified. Additional deviations from 
the s~^ decay of the cross section are due to logarithmic corrections which 
have their origin in the running coupling constant as and eventually in the 
evolution of the DAs (p^ (neglected in our calculation) . 



3.1 The K+-A Channel 

Figure R^ (left) shows the diquark-model predictions for the (scaled) photo 
production cross section s'^da/dt along with the few existing photoproductio, 
data at large- momentum transferlla and the outcome of the pure quark HSA 
(dash-dotted curve) . Whereas the DAs of proton and A have been kept fixed 
according to Eq. (0) we have varied the K+ DA. The solid and the dashed line 
represent results for the asymptotic (Eq. (0)) and the two- humped (Eq. (m) 
K+ DA, respectively, evaluated at pj^j^ = 6 GeV. The better performance 
of the asymptotic DA and the overshooting of the asymmetric DA is in line 
with the couiJusion drawn from the investigation of the pion-photon transition 
form factor L3 where, for the case of the pion, (j)^^ leads also to an overshoot- 
ing the data. Our findings have to be contrasted with those obtained within 
the pure quark-model calculation of photoproduction Q, where the asymptotic 
forms for both, baryon and meson DAs, give systematically larger results than 
the combination of very asymmetric DAs. However, the numerics of Rcf. 7 
must be taken with_some provisio. For Compton scattering off nucleons it has 
been demonstratedllj that the very crude treatment of propagator singularities 
adopted in Rcf. 7, namely keeping ie small but finite, may lead to deviations 
from the correct result which are as large as one order of magnitude. At this 
point we want to emphasize that we have paid special attention to a correct 
and numerically robust treatment of the propagator singularities which oc- 
cur in the convolution integral, Eq. (|l|). By carefully separating the singular 
contributions, exploiting delta functions, rewriting principal-value integrals as 
ordinary integrals plus analytically solvable principle-value integrals, we are 
able to do all the numerical integrations by means of very fast and stable 
fixed-point Gaussian quadratures. 

We have also examined the relative importance of various groups of Feyn- 
man graphs and found the 3-point contributions to be by far the most im- 
portant (cf. Fig. (left)). 4- and 5-point contributions amount to « 5% 
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Figure 3: Diquark-model predictions for photo- (left figure) and electroproduction (right 
figure) of the K+-A final state. Results are shown for fixed pZ. = 6 GeV. The experimental 
points are photoproduction data taken from Anderson et al.Hj. Left figure: differential cross 
section for 7p — » K+ A scaled by s^ vs coa/0cm); solid (dashed) line: dLquark-nmdel result 
for p and A DAs chosen according to Eq. (H), K"*" DA according to Eq. (ph (Eq. (M); dotted 
line; same as full line, tmt only three-point contributions taken into account; dash-dotted 
line; quark-model result El for the asymmetric p and A DAs of Ref. 13 and the two- humped 
K+ DA of Eq. Mj. Right figure: the spin-averaged cross section dau/dt (scaled by s^) vs 
cos(6cm) for transversely polarized photons with virtualities Q^ = (photoproduction limit), 
0.5, 1, 2, and 3 GeV^; dashes become shorter with increasing Q^. Proton and A DAs are 
chosen according to Eq. (0), the K+ DA according to Eq. (H). 



at 9cni — 90° and E^^^^ — 6 GeV as long as only da/dt is considered. Spin 
observables, on the other hand, are much more affected by 4- and 5-point con- 
tributions. A more detailed discussion of 7 p — > K"*" A (and also 7 p -^ K*^A) 
with a full account of calculational techniques and analytical expressions for 
the photoproduction amplitudes can be found in Ref. 16. 

The diquark-model predictions for a few electroproduction observables are 
depicted in Fig. ^ (right) and Fig. ^. In these plots we have concentrated on the 
asymptotic form, Eq. (^, of the K+ DA. Like in the case of photoproduction 
the large-s behaviour of the four cross section contributions is s^^ (modified 
by logarithmic factors) provided that Q"^ / s is kept fixed. For fixed photon 
virtuality Q^ and s — > 00, however, dai^/dt and dai/dt decay like s~^ and s~^, 
respectively. For fixed p^^y^ = 6 GeV the transverse cross-section contribution 
da\j/dt decreases with increasing Q^/s (cf. Fig. || (right)). The same holds for 
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Figure 4: Diquark-model predictions for electroproduction of the K"'"-A final state. 
Left figure: ratio of cross sections daj^/dau for longitudinally and transversely polarized 
photons, respectively, vs cos(9cm) for fixed p^^j^ = 6 GeV; dashes become shorter with in- 
creasing photon virtuality Q^ (= 0.5, 1, 2, 3 GeV'^). Right figure: dau/dt at fixed 9cm = 90° 
vs Q^/s for s = 10 (solid), 20 (long dashed), and 100 GeV^ (short dashed). 



dai^/dt (if (3^/s>0.04). At Q^/s = 0.04 doi^/dt is already larger than da^/dt 
(cf. Fig. El (left)). Since daj^/dt vanishes for Q^ ^- it thus has to rise very 
sharply at small values of Q'^/s. To get a better feeling for the Q^-dependence 
of dau/dt we have plotted this quantity for fixed scattering angle (0cm = 90°) 
and different values of s as function of Q^/s (cf. Fig. |j (right)). For exact 
s~^ scaling the curves for the three different values of s should coincide. The 
deviation from the s~^ scaling can mainly be ascribed to the running coupling 
constant Og. The dependence of dcr\j /dt (and likewise dai^/dt) on Q^/s can be 
roughly parameterized by means of a function ex (1 + cQ'^/s)^^ with the pa- 
rameter c between 2 and 2.5. The qualitative features, like scaling behaviour, 
Q'^/s dependence, and dominance of dui^/dt as compared to dau/dt remain, 
of course, unaltered if the asymptotic K+ DA is replaced by the asymmetric 
DA of Eq. (^. The angular dependence of cross-section ratios, like dai^/dajj, 
da^/dau, or dai/dij\j, however, exhibit a marked sensitivity on the choice of 
the K"*" DA. First experimental constraints on hard exclusive electroproduc- 
tion are to be expected from CEBAF and later on at even higher momentum 
transfers from ELFE (at DESY). 



10 



10 



^- 10^ 



> 




(IJ 




o 


10 


a 








1 . 




■n 




B 


10 



10'' 



10= 



; ♦ P\ab = 4GeV ; 


r A pT^,,b = 6Ge' 


1 


^ ^/^ 


" 



10' 



-0.8 



-0.4 0.4 
cos(e) 



10 



10 



> Q 

^ 7 



10" 



10= 



: ♦ pVb = 4GeV 


rA p\,b = 5GeV , 


• P\ab = 7-5GeV // 






: ^ 


y 


~ 1 



-0.8 



-0.4 0.4 
cos(0) 



0.8 



Figure 5: Diquark-model predictions for photoproduction of the K^-E" {left figure) and 
the 7r+-n {right figure) final states; solid (dashed) line: diquark-model results for baryon 
DAs chosen according to Eq. (Eh, K and n DAs according to Eq. (H) (Eq. (W))- Results are 
shown for fixed p^^^^ = 6 GeV. Scalar-diquark contributions are fully taken into account, 
four- and five-point contributions of vector diquarksn are neglected. Data are taken from 

Anderson et al.l'H. 



3.2 The K+-i;° and the 7r+-n Channel 

We have already mentioned that these channels differ from K+-A production 
as far as also V-diquarks are involved. Since the treatment of V diquarks is 
much more intricate than that of S diquarks we have, until now, only calculated 
the corresponding three-point contributions. From K+-A production and also 
from other applications of the diquark model we know, however, that spin- 
averaged cross sections are mainly determined by three-point contributions. 
Thus we do not expect a significant modification of our numerical results for 
the K^-S" and the tt^-u photoproduction cross sections if four- and five-point 
contributions of V diquarks are included. Like in the case of K+-A production 
we observe reasonable agreement with the K+-E° and the 7r+-n cross-section 
data if the asymptotic DA is taken for the K"*" and the 7r+, respectively (cf. 
Fig. H). The two- humped DAs, Eq. (Q), seem again to be in conflict with the 
data. Also for these two channels the pure quark HSA fails in reproducing the 
data. It is remarkable that the unpolarized differential cross sections for K+-A 
and K^-E" production are very similar in size, although the corresponding 
production mechanisms (via S and V diquarks, respectively) are quite different. 
We expect this difference to show up more clearly in spin observables. 
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4 Conclusions 

The predictions of the diquark model for 7P -^ K+ A, K+ S°, and 7r+ n look 
rather promising if the asymptotic form (ex x{l — x)) is taken for the me- 
son DAs. To the best of our knowledge the diquark model is, as yet, the 
only constituent scattering model which is able to account for the large-p_L 
photoproduction data. For electroproduction of the K+-A final state our re- 
sults are the first perturbative QCD predictions at all. With respect to future 
experiments it would, of course, be desirable to have more and better large 
momentum-transfer data on exclusive photo- and electroproduction. Polariza- 
tion measurements of the recoiling particle could help to decide, whether the 
perturbative regime has been reached already, or whether non-perturbative ef- 
fects (different from diquarks) are still at work. Spin observables, in general, 
could be very helpful to constrain the form of the hadron DAs. 
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